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Abstract. We provide a general and flexible approach to LIBOR mod- 
f*"^ ■ eling based on the class of affine factor processes. Our approach respects 

O^l ' the basic economic requirement that LIBOR rates are non-negative, and 

f*, . the basic requirement from mathematical finance that LIBOR rates 

' are analytically tractable martingales with respect to their own for- 

Lj^ ( ward measure. Additionally, and most importantly, our approach also 

leads to analytically tractable expressions of multi-LIBOR payoffs. This 
approach unifies therefore the advantages of well-known forward price 
models with those of classical LIBOR rate models. Several examples are 
added and prototypical volatility smiles are shown. We believe that the 
fV^ CIR-process based LIBOR model might be of particular interest for ap- 

0^ , plications, since closed form valuation formulas for caps and swaptions 

' ' are derived. 

1. Introduction 

Let To < . . . < Tj\r be the discrete tenor of maturity dates. LIBOR rates 
are calculated from the observable ratio of prices of zero-coupon bonds with 
. maturity T^,„i and = + 5 via 

si ^^(w- 1 

It is clear, due to the very nature of interbank loans, that LIBOR rates 
■ should be non- negative. Additionally, as a requirement from mathematical 

, finance, LIBOR rates should be martingales with respect to their own for- 

ward measure PT fc ; that is, when B(-,Tk) is considered as numeraire of the 

Bit T ) 

model, then discounted bond prices ( B(tTk) \ <t<T _ should be martin- 
gales. An additional basic requirement of the model is its tractability, since 
otherwise one cannot calibrate the model to the market data. Therefore the 
LIBOR rate processes L(-, T&) should have tractable stochastic dynam- 
ics with respect to their forward measure P^, for k = 1, . . . , N; for instance 
of exponential Levy type along the discrete tenor of dates Tq < . . . < T/v • 
Here the terminus "analytically tractable" is used in the sense that either 
the density of the stochastic factors driving the LIBOR rate process is known 
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explicitly, or its characteristic function. In both cases, the numerical evalu- 
ation, which is needed for calibration to the market, is easily done. 

In applications, the stochastic factors have to be evaluated with respect 
to different numeraires. In order to describe the dynamics with respect to 
a suitable martingale measure, for instance the terminal forward measure 
TPt N i we have to perform a change of measure. Usually this change of mea- 
sure destroys the tractable structure of L(-,Tk-\,Tk) with respect to its for- 
ward measure. This well-known phenomenon makes LIBOR market mod- 
els based on Brownian motions or Levy processes quite delicate to apply 
for multi-LIBOR-dependent payoffs: either one performs expensive Monte 
Carlo simulations or one has to approximate the equation (the keyword here 
is "freezing the drift", see Siopacha and Teichmann 2010[ ). 



In order to overcome this natural intractability, forward price models 
have been considered, where the tractability with respect to other forward 
measures is pertained when changing the measure. Hence, modeling forward 
prices F(-,Tk-i,Tk) = 1 + 5L(-,Tk) produces a very tractable model class, 
however, negative LIBOR rates can occur with positive probability, which 
contradicts any economic intuition. 

In this work we propose a new approach to modeling LIBOR rates based 
on affine processes. The approach follows the footsteps of the forward price 
model, however, we are able to circumvent the drawback of traditional for- 
ward price models: in our approach LIBOR rates are almost surely non- 
negative. Furthermore, closed form valuation formulas - using Fourier trans- 
forms - for caplets and swaptions, or any other multi-LIBOR payoff can be 
derived, hence the calibration and evaluation of those models is fairly simple. 
In fact, the model remains analytically tractable with respect to all possible 
forward measures which can be chosen. 

A particular feature of our approach is that the factor process is a time- 
homogenous Markov process when we consider the model with respect to 
the terminal measure ~R?t n - With respect to forward measures the factor 
processes will show time-inhomogeneities due to the nature of the change of 
measure. When we compare our approach to an affine factor setting within 
the HJM-methodology, we observe that in both cases one can choose - with 
respect to the spot measure in the HJM setting or with respect to the termi- 
nal measure in our setting - a time-homogeneous factorprocess. LIBOR-rates 
have in both cases a typical dependence on time-to-maturity T/v — t. 

The remainder of the article is organized as follows: in Section [2] we for- 
mulate basic axioms for LIBOR market models. In Section[3]we recapitulate 
the literature on LIBOR models. In Section 0] we introduce affine processes 
which are applied in Section [5] for the construction for certain martingales. 
In Section [6] we present our new approach to LIBOR market models, which 
is applied in Section [7] to derivative pricing. In Section [8] several examples, 
including the CIR-based model, are presented and in Section [9] we show 
prototypical volatility surfaces generated by the models. 



2. Axioms 



Let us denote by L(t, T) the time-i forward LIBOR rate that is settled at 
time T and received at time T + 5; here T denotes some finite time horizon. 
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The LIBOR rate is related to the prices of zero coupon bonds, denoted by 
B{t,T), and the forward price, denoted by F(t,T,T + 6), by the following 
equations: 

1 + SL(t, T) = B B ^ T T } 8) = F(t, T,T + S). (2.1) 

One postulates that the LIBOR rate should satisfy the following axioms, 
motivated by economic theory and arbitrage pricing theory. 

Axiom 1. The LIBOR rate should be non-negative, i.e. L(t,T) > for all 
< t < T. 

Axiom 2. The LIBOR rate processes should be a martingale under the cor- 
responding forward measure, i.e. L(-,T) G A4(Pt + $). Moreover, the LIBOR 
rate processes should be analytically tractable with respect to as many for- 
ward measures as possible; minimally, closed-form valuation formulas should 
be available for the most liquid interest rate derivatives, i.e. caps and swap- 
tions, so that the model can be calibrated to market data in reasonable time. 

Furthermore we wish to have rich structural properties: that is, the model 
should be able to reproduce the observed phenomena in interest rate mar- 
kets, e.g. the shape of the implied volatility surface in cap markets or the 
implied correlation structure in swaption markets. 

Here, a "closed-form" valuation formula refers to the explicit knowledge of 
the density or the explicit knowledge of the characteristic function, whence 
Fourier-based methods for the pricing of European options apply. 

We refer to the first axiom as an economic requirement deduced from 
the structure of interbank markets. The second axiom is referred to as a 
requirement from mathematical finance. 

3. Existing approaches 

There are several approaches developed in the literature attempting to 
fulfill these axioms and practical requirements. We briefly describe below 
the two main approaches, namely the LIBOR market model and the for- 
ward price model, and comment on their ability to fulfill these axioms and 
requirements. 

Approach 1. In LIBOR market models, developed in a series of articles 



by Sandmann et al. (1995) Miltersen et al. (1997) Brace et al. (1997) and 



Jamshidian (1997) , each forward LIBOR rate is modeled as an exponential 
Brownian motion under its corresponding forward measure; this model pro- 
vides a theoretical justification for the common market practice of pricing 
caplets according to Black's futures formula (jBlack 1976|) . i.e. assuming that 
the forward LIBOR rate is log-normally distributed. Several extensions of 
this framework have been proposed in the literature, using jump-diffusions, 
Levy processes or general semimartingales as the dr iving motion (cf. e.g. 



IGlasserman and Kou 2003], lEberlein and Ozkan 2005| IJamshidian 1999]) , or 



incorporating stochastic volatility effects (cf. e.g. Andersen and Brotherton- 
Ratcliffe MM- 

We can generically describe LIBOR market models as follows: on a sto- 
chastic basis consider a discrete tenor of dates (7fc)o<fc<jv, forward measures 
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JPx k associated to each date and appropriate volatility functions A(-, T k ) also 
associated to each date; let H be a semimartingale starting from zero, with 
predictable characteristics (B,C,v) or local characteristics (b,c,F) under 
the terminal measure TPt n , driving all LIBOR rates. Then, the dynamics of 
the forward LIBOR rate with maturity T k is 

L(t,T k ) = L(0,T fc )exp (j b(s,T k )ds + J X(s,T k )dHj k+1 \ > 0, (3.1) 

where H Tk+1 denotes the martingale part of the semimartingale H under 
the measure TPx k+1 , and the drift term is 

b(s,T k ) = -^X(s,T k ) 2 c s 

- J (e x ^> - 1 - X(s,T k )x)Fl k+1 (dx), (3.2) 

k 

ensuring that L(-,T k ) 6 M(TPx k+1 )- The semimartingale H Tk+1 has the 
IFVfc+i -canonical decomposition 

t t 

H? k+1 = j yfedW? k+1 + j j x([i H ~ v Tk+1 )(ds,dx), (3.3) 

OR 

where the Pt^ -Brownian motion is 



* ' N 



1 e /ff;: T l ^,T ; ))^d g , (3.4) 

\l=k+l 







and the Pr fe+1 -compensator of \x H is 



(3.5) 



As an example, the classical log-normal LIBOR model is described in this 
context by setting (6, c, F) = (0, cr 2 , 0). 

Now, let us discuss some consequences of this modeling approach; clearly 
H remains a semimartingale under any forward measure, since the class of 
semimartingales is closed under equivalent measure changes. However, any 
additional structure that we impose on the process H to make the model 
analytically tractable will be destroyed by the measure changes from the 
terminal to the forward measures, as the random terms -, tttt?-^ \ enter- 
ing into eqs. (|3.4|) and (|3.5|) clearly indicate. For example, if H is a Levy 
process under TPt n , then H Tk+1 is not a Levy process (not even a process 
with independent increments) under Py H1 . Hence, we have the following 
consequences: 

(1) if H is a continuous semimartingale, then caplets can be priced in 
closed form, but not swaptions or other multi-LIBOR derivatives; 
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(2) if H is a general semimartingale, then even caplets cannot be priced 
in closed form. 

Moreover, the Monte Carlo simulation of LIBOR rates in this model is 
computationally very expensive, due to the complexity evident in eqs. (|3.4p 
and (|3.5p . Expressing the dynamics of the LIBOR rate in (|3.1|) under the 
terminal measure leads to a random drift term, hence we need to simulate 
the whole path and not just the terminal random variable. More severely, 
the random drift term of e.g. L(-,Tk) depends on all subsequent LIBOR 
rates L(-,T[), k + 1 < I < N; hence, to simulate L(-,Tfe) we must previously 
simulate - or restore from the memory - the paths of all subsequent LIBOR 
rates L(-, 7]) for all k + 1 < I < N. Indeed, the dependence of each rate on all 
rates with later maturity can be represented as a strictly (lower) triangular 
matrix. 

Of course, some remedies for the analytical intractability of the LIBOR 
market model have been proposed in the literature. The common practice is 
to replace the random terms in (|3,4|) and (|3.5|) by their deterministic initial 
values, i.e. to approximate 

SiL(t-,Ti) _ SiL(0,Ti) . 



1 + W-T,) 1 + ^(0,1})' 

this is usually called the "frozen drift" approximation. As a consequence, the 
structure of the process H will be - loosely speaking - preserved under the 
successive measure changes; for example, if H is a Levy process, then H Tk+1 
will become a time-inhomogeneous Levy process (due to the time-dependent 
volatility function). Hence, caps and swaptions can be priced in closed form. 
However, this is an "ad hoc" approximation, with no theoretical justification, 
and error estimates are not available. Moreover, simulation results show that 
the quality of the approximation deteriorates for more complex multi-LIBOR 
products. 



More recently, Siopacha and Teichmann (2010) and Papapantoleon and 



Siopacha (2 009p have developed Taylor approximation schemes for the ran- 
dom terms entering (|3.4|) and (|3.5p using perturbation-based techniques. 
This method offers approximations that are much more precise than the 
"frozen drift" approximation (|3,6p . while at the same time being significantly 
faster than simulating the actual dynamics; moreover, it offers a theoretical 
justification for the "frozen drift" approximation as the zero-order Taylor ex- 
pansion. However, this method is based on Monte Carlo simulations, hence 
is not fast enough for calibration in real time. 

Therefore, this approach satisfies Axiom 1. As far as Axiom 2 is con- 
cerned, LIBOR rates are analytically tractable only under their own for- 
ward measure and only if the driving process is continuous; LIBOR rates 
are not tractable with respect to any other forward measure. Therefore caps 
can (possibly) be priced in closed form, but not swaptions or more compli- 
cated multi-LIBOR derivatives. 

Remark 3.1. Additionally it is econometrically not desirable to model LI- 
BOR rates as exponentials of processes with independent increments. How- 
ever, we admit that this is a minor point. 
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Approach 2. In the forward price model proposed by Eberlein and Ozkan 



and Kluge (2005) , the forward price - instead of the LIBOR rate 
- is modeled as an exponential Levy process or semimartingale. Consider 
a setting similar to the previous approach: (Tfc) <fc<jv is a discrete tenor 
of dates, are forward measures and H denotes a semimartingale with 
characteristics (B,C,v) under the terminal measure TPt n , where Hq = 0. 
Then, the dynamics of the forward price F(-, T k , T k+ i), or equivalently of 
1 + SL(;T k ), is given by 

(t t > 

J b(s,T k )ds + J X(s,T k )dHl k+1 
/ 

(3.7) 

where H Tk+1 denotes the martingale part of H under the measure PT fc+1 
and the drift term b(s,T k ) is analogous to (|3.2p . ensuring that L(-,T k ) G 
A4(TPx k+1 )- The semimartingale H Tk + 1 has the PT fc+1 -canonical decomposi- 
tion 

t t 
Hj k+1 = J ^c~ s dwl k+1 + J J x((i H -v T ^)(ds,dx), (3.8) 

R 

where the Pj- fc+1 -Brownian motion is 



lF, r ' +1 =W t - I I > ' A(t,T,) I ^ds, (3.9) 



/ f E A(*,T,)) 

J Q \l = k + l / 



and the PT fc+1 -compensator of fi H is 

u Tk+1 (ds , dx) = exp I x A(s,T/) I v(ds,dx). (3.10) 




Now, we can immediately deduce from (|3.9p and ()3.10p that the structure 
of the process H under Ptv is preserved under any forward measure Pr fc+1 . 
For example, if H is a Levy process under Ptjv then it becomes a time- 
inhomogeneous Levy process under , if the volatility function is time- 

dependent; we can also deduce that the measure change from the terminal to 
any forward measure is an Esscher transformation (cf. Kallsen and Shiryaev 



As a result, the model is analytically tractable, and caps and swaptions 
can be priced in closed form (similarly to an HJM model). However, negative 
LIBOR rates can occur in this model, since forward prices in those models 
are usually positive but not necessarily greater than one. Hence Axiom 1 
is violated, but Axiom 2 is satisfied in the best possible way, since LIBOR 
rate processes are analytically tractable with respect to all possible forward 
measures. 

Remark 3.2. The forward price model can be embedded in the HJM frame- 



work with a deterministic volatility structure; cf. Kluge (2005, §3.1.1. 
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The two modeling approaches we have just reviewed might appear similar 
at first sight, but they actually differ in quite fundamental ways - apart from 
the considerations regarding Axioms 1 and 2. 

On the one hand, the distributional properties are markedly different; in 
the LIBOR market model - driven by Brownian motion - LIBOR rates are 
log-normally distributed, while in the forward price model - again driven by 
Brownian motion - LIBOR rates are, approximately, normally distributed. 
Although there seems to be no consensus among market participants on 
which assumption is better, it is worth pointing out that in the CEV model 
- where for /3 — > the law is normal and for (3 — > 1 the law is log-normal - 
a typical value for market data is /3 ~ 0.4. 

On the other hand, changes in the driving process affect LIBOR rates in 
the LIBOR model and the forward price model in a very different way; see 



also Kluge (2005, pp. 60) Assume that in a small time interval of length dt 
the driving process changes its value by a small amount A. Then, in the 
LIBOR model we get: 

Lit + dt, T) m L(t, T) + A ■ L(t, T) + C(A 2 ), (3.11) 

while in the forward price model we get: 

L{t + dt, T) « L(t, T) + ^ + A • L{t, T) + 0{A 2 ). (3.12) 

o 

Hence, in the LIBOR model changes in the driving process affect the rate 
roughly proportional to the current level of the LIBOR rate; in the forward 
price model, changes do not depend on the actual level of the LIBOR rate. 

Aim: we would like to construct a "forward price" -type model with positive 
LIBOR rates; that is, we want an model that respects Axioms 1 and 2. 

A first idea would be to search for a process that makes the martingale 
in ()3.7p greater than one, hence guaranteeing that LIBOR rates are always 
positive. However, such an attempt is doomed to fail since one has to demand 
in equation (j3.7|) that 



j b(s,T k )ds + J \(s,T k )dHj k+1 >0 



with respect to the foward measure (or any other equivalent measure). This 
reduces the class of available semimartingales considerably and restricts the 
applicability of the models. We show in Section [5] an alternative construction 
with rich stochastic structure. 



4. Affine processes 

Let (O, F, F, P) denote a complete stochastic basis, where F = {Tt)t€[o,T\ > 
and let < T < oo denote some, possibly infinite, time horizon. We consider 
a process X of the following type: 

Assumption (A). Let X = (Xt)o<t<T be a conservative, time-homogene- 
ous, stochastically continuous Markov process taking values in D = M>q, 
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and (JP x )xeD a family of probability measures on (f2, F), such that Xq = x 
P^-almost surely, for every x £ D. Setting 

1 T : = | it £ M d : E x [e (u ' X:r) ] < oo, for all x £ I)| , (4.1) 

we assume that 
(i) £ Z° ; 

(ii) the conditional moment generating function of Xt under TP X has 
exponentially-affine dependence on x; that is, there exist functions 
4>t\u) : [0, T]xI T ^K and ip t (u) : [0, T]xI T ^ M d such that 

E x [exp(u,X t )] =exp(^(u) + (^(u),x)), (4.2) 

for all (t,u,x) £ [0,T] xl T x D. 

Here "•" or (•, •) denote the inner product on W 1 , and E x the expectation 
with respect to P x . 

Stochastic processes on ]R^ with the "affine property" f|4.2j) have been 
studied since the seventies as the continuous-time limits of Galton- Watson 



branching processes with immigration, cf. Kawazu and Watanabe (1971) 
More recently, such processes on the more general state space M™ x W 1 
have been studied comprehensively, and with a view towards applications in 
finance, by Duffie, Filipovic and Schachermayer (|2003p . We will largely follow 



their approach, complemented by some results from Keller-Ressel (2008) 



By Theorem 3.18 in Keller-Ressel (2008) , the right hand derivatives 



| d 

F{u) := —\ t=Q+ (t) t {u) and R{u) := —\ t=0+ ?p t (u) (4.3) 
exist for all u £ It and are continuous in u, such that X is a 'regular affine 



process' in the sense of Duffie et al. (2003) Moreover, F and R satisfy Levy 



Khintchine-type equations; it holds that 

F(u) = (b,u) + J (e^' u > - l))m(d£) (4.4) 



D 

and 



Ri(u) = (fr,u) + + / {^' u) ~ 1 - (u,h\t)))^(dO, (4.5) 



D 

where (b, m, cti, fy, ^i)i<i<d are admissible parameters, and h l : R^, — > M. d 
are suitable truncation functions, defined coordinate- wise by 

hUO--=\ ) £ s k u t % - foralU£M| ,i£{l,...,4, (4.6) 

I^Xisfc)) K — l 

with x( z ) an y bounded Borel function that behaves like z in a neighborhood 
of 0, such as or zl{| z |< 1} . 

The parameters (b,m,ati, Pi, Hi)\<i<d have the following form: (Pi)i<i<d 
and b are Revalued vectors, (ai)i<j<d are positive semidefinite real d x d 
matrices, and m and (fJ*i)i<i<d are Levy measures on R> , satisfying addi- 
tional admissibility conditions; writing I = {1,. . . ,d}, these conditions are 
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riven, according to Duffie et al. (2003) , by 

b g M| (4.7) 

ft(fc) e^oVfce /\W and G E (4.8) 

a i(fe0 =0 if fc € or I G (4.9) 



m({0}) = and y (|f | A l)m(d£) < oo (4.10) 
D 

where |£| = J2i \d\ for £ G R d ; and, for alH G / 

Mi({0})=0 and /" [(|£ A{i} | + |&| 2 ) A l] Md£) < oo. (4.11) 
D 

The time-homogeneous Markov property of X implies the following con- 
ditional version of (14.21): 



TE x [exp(u,X t+s )\F s ] =exp((f> t (u) + (Mu),X s )), (4.12) 

for all < t + s < T and u G It- Applying this equation iteratively, it is 
seen that the functions (ft and tft satisfy the semi-flow property 

<ftt+s{u) = (ftt{u) + 4> s {tft t {u)) 
lftt+s(u) = lft s (ll>t(u)) 

for all < t + s < T and u G It, with initial condition 

(fto(u) = and ifto(u) = u; (4-14) 



see also Lemma 3.1 in Duffie et al. (2003) and Proposition 1.3 in Keller- 
Ressel (|2008l) . 

Differentiating the flow equations (and using the existence of (|4.3[) ) we 
arrive at the following ODEs (the generalized Riccati equations) satisfied by 
(ftt and ift t - 

= F(Mu)), 0o («) = 0, (4.15a) 
^(u) = R(M»)), Mu) = u, (4.15b) 



for (t,u) G [0,T]xI T ; cf. Duffie et al. (2003, Theorem 2.7) If ift t (u) stays in 



I T for all t G [0,T], it is a unique solution. Note that if the jump measures 
m and \x are zero, F(u) and each Ri(u) are quadratic polynomials, whence 
the differential equations degenerate into classical Riccati equations. 

Finally, let us mention, that any choice of admissible parameters satisfying 
(|4.7p - (|4.11|) . and corresponding functions F and R, gives rise to a uniquely 
defined affine process, whose moment generating function can be calculated 
through the generalized Riccati equations (|4.15|) . 

Remark 4.1. We mention here the following examples of one-dimensional 
processes, satisfying Assumption (A): 

(1) Every Levy subordinator with cumulant generating function k(u) 

and finite exponential moment; it is characterized by the functions 

F(u) = k(u) and R(u) = 0. 
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(2) Every OU-type process (cf. ISato 19991 section 17) driven by a Levy 
subordinator with finite exponential moment; such a process is char- 
acterized by F(u) = k(u) and R(u) = (3u, with ft G R. 

(3) The squared Bessel process of dimension a (cf. IRevuz and Yor 1999|, 
Ch. XI), characterized by F{u) = au and R(u) = 2u 2 , with a > 0. 

Finally, we will later need the following results; let us denote by (e ; ),<d 
the unit vectors in M. d and let inequalities involving vectors be interpreted 
component- wise. 

Lemma 4.2. The functions 4> and i/j satisfy the following: 

(1) &(0) = Vt(0) = for all t G [0, T]. 

(2) It is a convex set; moreover, for each t G [0, T], the functions It B 
u i—)- 4>t{u) and It 3 u \-t <frt(u) are (componentwise) convex. 

(3) </>t(-) and ipt(') are order-preserving: let (t,u),(t,v) G [0,T] x It, 
with u < v. Then 

Mu) < <k(v) and %Ih{u) < M^)- (4-16) 

(4) ipt(') is strictly order-preserving: let (t,u), (t,v) G [0,T] xI T , with 
u < v. Then ipt{u) < ipt(v). 

Proof From f|4.4|) and (|4.5p it is immediately seen that F(0) = R(0) = 0. 
Thus (j>t(0) = V'i(O) = are solutions to the corresponding generalized Ric- 
cati equations (|4.15p . Moreover, G I T , such that the solutions are unique, 
showing claim (1). Let u, v G W 1 and A G [0, 1]. By Holder's inequality 

TE x [exp((\u + (l-\)v,X t ))] <TE x [e^] X -TE x [e^f- X) , (4.17) 

where both sides may take the value +oo. Taking logarithms on both sides 
shows that, for all t G [0, T], (fit(-) and V't(') are (componentwise) convex 
functions on M d , taking values in the extended real numbers 1R U {+00} . 
This implies in particular that It is convex, and that the restrictions of 
4>t(-) and ipt(') toXf are finite convex functions, showing claim (2). Following 



Keller- Ressel (2008, Proposition 1.3(vii)) , we have that for u < v 



E^e^)] <E* !><"•*>] <oo, 

for all x G K> . Now, using the affine property of the moment generating 
function we get 

<Pt(u) + < <j)t{v) + (Mv),x), (4.18) 

whereby inserting first x = and then x = Ce,, for C > arbitrarily 
large, yields claim (3). Consider the Riccati differential equation (|4.15b|) . 
satisfied by ipf By Keller-Ressel (2008) , Lemma 4.6, R(u) is quasi-monotone 
increasing; moreover, it is locally Lipschitz inZ£. A comparison principle for 
quasi-monotone ODEs (cf. lWalter 19961 Section 10. XII) yields then directly 
that u < v implies tpt(u) < ipt(v) for all t G [0,T]. □ 

The above results on affine processes can be extended to the case when 
the time-homogeneity assumption on the Markov process X = (Xt)o<t<T is 
dropped, see Filipovic (2005) The conditional moment generating function 
then takes the form 

E x [exp(u,X r }\ T s ] = exp (<f> 8 , r (u) + (ip s , r (u),X s )), (4.19) 
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for all (s, r, u) such that < s < r < T and u G It, with (j) Str (u) and ip sr (u) 
now depending on both s and r. Assuming that X satisfies the 'strong 
regularity condition' (cf. Filipovic 2005, Definition 2.9), <f> s , r (u) and ip sr (u) 
satisfy generalized Riccati equations with time-dependent right-hand sides: 

d 

.(u) = F(s,tp St r(u)), 



ds 



R(s,lp S)r (u)), 



4>r,r(u) 



0. 



it, 



(4.20) 
(4.21) 



for all < s < r < T and u G It- 



5. Constructing Martingales ^ 1 

In this section we construct martingales that stay greater than one for all 
times, up to a bounded time horizon T, that is, from now on < T < oo. 
The construction is a "backward" one, and utilizes the Markov property of 
affine processes. 

Theorem 5.1. Let X be an affine process satisfying Assumption (A), and 
let u G It- The process M" = (M")o<t<T defined by 

= exp {<h-M + (^ T - t (u),X t )), (5.1) 

is a martingale. Moreover, ifu G X^nR^Q then M" > 1 a.s. for allt G [0, T], 
for any X G M| . 

Proof. First, we show that M u is a martingale; for all u G It holds that 

E X [M%] =E x .[e< u ' Xr >] <oo. 

Moreover, using (|4.14p and (|4.12p . we have that: 

JE x [M^\T t ] =M x [esxp{u t X T )\J r t] 

= exp {<t> T -t{u) + (iJT-t(u),X t )) = M". 

Regarding the assertion that Mf > 1 for all t G [0, T], it suffices to note 
that if u G It H B^o> t^ 1611 is the conditional expectation of a process 
greater than, or equal to, one, i.e. 

M" = E x [exp(n, A T )|Ji] > 1, (5.2) 

hence greater than, or equal to, one itself. □ 

Remark 5.2. Actually, the same construction would create martingales for 
any Markov process X on a general state space, cf. Appendix [XJ However, 
taking the positive orthant as state space guarantees that the martingales 
stay greater than one; moreover, taking an affine process as the driving 
motion provides the appropriate trade-off between rich structural properties 
and analytical tractability. 

Remark 5.3 (Levy processes). Assume that the affine process X is actually 
a Levy subordinator, with cumulant generating function k. Then, we know 
that 



4>t(u) = t ■ k(u) and tpt{u) = u. 



(5.3) 
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Hence, the exponential martingale in (|5.ip takes the form: 

= exp (<j>T-t(u) + ^T-t{u),X t )) 

= exp((T-t)K(u) + (u,X t )), (5.4) 

which is a martingale by standard results for Levy processes. Moreover, for 
u G Zt, since k : It B^o an d T — t > 0, we get that M" > 1 for all 
t€[0,T\. 

Remark 5.4. On the other hand, this calculation shows that this model 



will contain the Levy forward price model of Eberlein and Ozkan (2005) and 



Kluge (2005) as a special case, if we consider a time-inhomogeneous affine 
process with state space M. d as driving motion. Of course, in that case the 
martingales M u will not be greater than one. 

Note that there is still some ambiguity lurking in the specification of the 
martingale M u : consider a d-dimensional driving process X, from which the 
martingale M u is constructed; let c be a positive semidefinite d x d matrix, 
and d its transpose. Define X = c-X, and let M" be the corresponding mar- 
tingale. It is easy to check that if X is an affine process satisfying condition 
A, then so is X. It holds that 

Mf u = [ exp(c'u, X T ) \T t ] = E x [ exp(u, cX T )\F t ] = M?, 

showing that in terms of the martingales M", a (positive) linear transfor- 
mation c of the underlying process X is simply equivalent to the transposed 
linear transformation d of the parameter u. In order to avoid this ambigu- 
ity in the specification of the martingale M u , we will fix from now on the 
initial value of the process X at some strictly positive, canonical value, e.g. 
1 = (1,...,1). 

Finally, the following definition will be needed later. 

Definition 5.5. For any process X = (X t )o<t<T satisfying Assumption 
(A), define 

7x = sup Ei[e (lt ' XT) ]. (5.5) 



6. LIBOR MODELING: A NEW APPROACH 

Now, we describe our proposed approach to modeling LIBOR rates that 
aims at combining the advantages of both the LIBOR and the forward price 
approach; that is, a framework that produces non-negative LIBOR rates in 
an analytically tractable model. 

Consider a discrete tenor = To < T\ < Ti < ■ ■ ■ < T/v = T and an initial 
tenor structure of non-negative LIBOR rates L(0, T/%), k £ {1, . . . ,N}. We 
have that discounted traded assets (bonds) are martingales with respect to 
the terminal martingale measure, i.e. 

frS e M(W Tn ), for all k 6 {1, . . . , N - 1}. (6.1) 
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Therefore, our ansatz is to model quotients of bond prices using the martin- 
gales M u as follows: 

f, ( *'lH = (6.2) 
B(t,T N ) 1 

= M; n ~\ (6.3) 



jB(t,7V_i) _ 



B{t,T N ) 

for all t 6 [0, Ti], . . . ,i G [0, T/v-i] respectively. Hence, the initial values of 
the martingales M Uk must satisfy: 

= exp (Mu k ) + <V*(«fc).*» = (6.4) 

for all k G {1, . . . ,N — 1}. Obviously we set itjv = 0« M U " = = 1. 

Next, we show that under mild conditions on the underlying process X, 
an affine LIBOR model can fit any given term structure of initial LIBOR 
rates through the parameters u±, . . . , lijv- 

Proposition 6.1. Suppose that L(0, Ti), . . . , L(0, T/v) is a tenor structure 
of non-negative initial LIBOR rates, and let X be a process satisfying as- 
sumption (A), starting at the canonical value 1. T/ie following hold: 

(1) If jx > B(0,Ti)/B(Q,Tn), then there exists a decreasing sequence 
ui > U2 > ■ ■ • > uat = in Tt n M>0) smc/i i/tai 

= £(0,?fc) y aHJfe€{l,...,JV}. (6.5) 

In particular, if jx = oo, then the affine LIBOR model can fit any 
term structure of non-negative initial LIBOR rates. 

(2) If X is one- dimensional, the sequence {uk)k£{i....,N} * s unique. 

(3) If all initial LIBOR rates are positive, the sequence (v>k)ke{i,...,N} ^ s 
strictly decreasing. 

Proof. The non-negativity of (initial) LIBOR rates clearly implies that 

fl(0,Tj) > B(0,T 2 ) > > B(0,T N ) _ 
5(0, T N ) ~ B(0,T N ) B(0, TV) 

Moreover, if the initial LIBOR rates are positive the above inequalities be- 
come strict. Now let e > 0, small enough such that 7x — e > b }q' jT J\ ■ Clearly, 
by the definition of 7x , we can find some u + > such that 

B(0,Ti) 



Ei 

Define now 



> lx - e > 



B(0, T N ) • 



/ : [0, 1] ->■ R >0 , ^Ej e K "+*) =M% U+ . (6.6) 

By monotone convergence, / is an increasing function; in addition Fatou's 
Lemma shows that / is lower semi-continuous. But any increasing, lower 
semi-continuous function is actually continuous. Thus / is a continuous, 
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increasing function satisfying /(0) = 1 and /(l) > b(0 Tn) ■ Consequently, 
there exist numbers = £n < 6v-i < • • • < £i < 1, such that 

/(&) = Mt + = ||^, foralUG{l,...,iV}. 

Setting u k = CkU+, we have shown ()6.5j) . By Lemma l4.2| /(£) is in fact 
a strictly increasing function. If also the (quotients of) bond prices satisfy 
strict inequalities, we deduce that the sequence (w&)fce{i,...,.iV} is strictly de- 
creasing, showing claim (3). Finally, if X is one-dimensional, then 1^ (1 M^o 
is just a sub-interval of the positive half-line; thus any choice of u+, will lead 
to the same parameters u k , showing (2). □ 

Forward prices have the following dynamics: 
B(t,T k ) B(t,T k ) B(t,T N ) Mt k 



B(t,T k+1 ) B(t,T N ) B(t,T k+1 ) M? k+1 
= exp \ 6 TN - t (u k ) - <j)T N -t(u k+1 ) 

+ (lpT N -t(u k ) - ^T N _t{u k+ i),X t ) 

= exp (A TN -t(u k ,u k+1 ) + (B TN - t (u k ,u k+1 ),X t )), (6.7) 

where we have defined 

Ar N -t(uk,Uk+i) ■= <htf-t( u k) ~ <h N -t( u k+i), (6.8) 
B TN _ t {uk, u k+1 ) := ipT N -t(u k ) - ip TN _ t (u k+1 ). (6.9) 

Using Proposition I6.1l fl) and Lemma I4.2f 3). we immediately deduce the 
following result, which shows that Axiom 1 is satisfied: 

Proposition 6.2. Suppose that L(0, T\), . . . , L(0, T/v) is a tenor structure 
of non-negative initial LIBOR rates, and let X be a process satisfying as- 
sumption (A) . Let the bond prices be given by (|6.2p - (|6.3|) and satisfy the 
initial conditions (|6.4|) . Then the LIBOR rates L(t,T k ) are non-negative 
a.s., for all t £ [0, T k ] and k G {1, . . . , N}. 

Moreover, forward prices should be martingales with respect to their cor- 
responding forward measures, that is 

3' : J k \ G MQP Tk+1 ), for alU € {1, . . . ,N - 1}; (6.10) 
B{-,l k+ i) 

this we can easily deduce in our modelling framework. Forward measures are 
related to each other via forward processes, hence in the present framework 
forward measures are related to one another via quotients of the martingales 
M u . Indeed, we have that 



dP Tfc 



dP Tfc 



+i 



F{t,T k ,T k+1 ) = B{0,T k+1 ) _M, 
T t F(0,T k ,T k+1 ) B(0,T k ) X M i 



(6.11) 



for any k E {1, . . . ,N — l},t £ [0,T k ]. Then, using Proposition III. 3. 8 in 



Jacod and Shiryaev (2003) we can easily deduce that L(-,T k ) is a martingale 
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under the forward measure lPj- k+1 , since the successive densities from ET fc+1 
to yield a "telescoping" product and a Py M martingale. We have that 



l + SL(;T k ) 



B(;T k ) 
B(;T> 



since 



N-l 

t n 

l=k+l 



fc+ly 



M u k gM (jp Tn ) 



(6.12) 



(6.13) 



by the construction of the model. 

In addition, we get that the density between the Py fe -forward measure 
and the terminal forward measure P^ is given by the martingale M Uk , as 
the defining equations (|6.2|) - (|6.3p already dictate; we have 



B(0,T N ) 

T t B(0,T k ) " B(t,T N ) B(0,T k ) 



B(0,T N ) B(t,T k ) 

X 



x Mf k 



(6.14) 



this we can also deduce by expanding the densities between P T k and P^ . 

Next, we wish to show that the model structure is preserved under any 
forward measure as well; indeed, changing from the terminal to the forward 
measure, X becomes a time-inhomogeneous Markov process, but the affine 
property of its moment generating function is preserved. This means that 
Axiom 2 is satisfied in full strength: X will be a time-inhomogeneous affine 
process under any forward measure; to show this, we calculate the condi- 
tional moment generating function of X r under the forward measure TPr k , 
and get that 



Ep 7 



M. 



1 



M. 
exp 



«rE P% [exp ((f) TN _ r (u k ) + (ip TN _ r (u k ),X r ) + (v,X r ))\T s 

(U k ) - {tpT N -s(Uk),X s ) + (j)T N -r(u k )) 



l N -s\ 



x Ep Tjv [exp ((ip TN - r (u k ) + v,X r )) \T s ] 

= exp (^> TN _ r (u k ) - 4>T N -s{uk) + <j)r-s(^T N -r(u k ) + v) 
+ {lp r - s (lpT N -r(u k ) + v) - ^ TN _ s (u k ), X s )j 
= exp ($ r -. s (lpr N - r (u k ) + V) - <j) r - s (lpT N -r(Uk)) 

+ (l/j r - s (lpT N -r(u k ) + v) - 1p r - s (lpT N -r(u k )),X s 



(6.15) 



which yields the affine property of X under the forward measure PT fc > for 
any k € {1, . . . , JV — 1}. In particular, setting s = 0, r = t, we get that 



E Prt [e^ Xt) ] = exp U{v) + x> 



(6.16) 
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where 

4{v) := cMV'twK) + v) - (f>t(i/JT N -t(u k )), (6.17) 

i/j%(v) := M^T N -t(u k ) +v)- M^T N -t(u k )), (6.18) 

showing clearly that the measure change from TPr k to TPt n is an expo- 
nential tilting (or Esscher transformation). Furthermore, we can calculate 
from (|6.15p the functions F k (r,v) and R k (r,v), characterizing the time- 
inhomogeneous affine process X under the forward measure IPr fe : 

d d 

F k (r,v) = - — (p r - s {i)T- T {uk) +V)\ s=r + —<j)r-s{^T~r{u k ))\ s=r 

= F {^ T -r{u k ) +V)-F (^T-rK)) , (6.19) 

and 

d d 

R k (r,v) = -— 1p r - s (lpT-r(u k ) +v)\ s=r + —*p r - s (lp T -r(u k ))\ s=r 

= R (iJT-r{Uk) +V)-R {^ T -r{u k )) . (6.20) 

Note that the moment generating function in (|6.16p is well defined for all 
ti£l l , where 

Z k = {v G R d : v + iP TN -t{u k ) el T ,te [0, T k }}. 

Moreover, we would like to calculate the moment generating function for 
the dynamics of forward prices under their corresponding forward measures. 
Let us use the following shorthand notation for 

B(t,T k ) _ ^A k +B k -X 



B(t,T k+ i 

where 

A k := A TN -t(u k ,u k+ i) = (f>T N -t(u k ) - 4>T N -t(u k+ i), 
B k := B TN -t(uk,u k+1 ) = ipT N -t{u k ) - ipT N -t{uk+i)- 
for any k G {1, . . . , N — 1}. Then, using (|6.16p we get that 
Ep T [e 

T k+1 L 



e**-^*^*, (6.21) 

(6.22) 



jv(A k +B k -X t )] 

k + 1 



; Ep t [. 

T k+1 L 

= e vAk exp U>t(ipT N -t(uk+i) + vB k ) - (j} t (tpT N -t(u k +i)) 

+ {A(^T N -t(u k +i) + vB k ) - ip t (ipT N - t (u k+ i)),x 
B(0,T N ) / , , . , . , . , 

= r X exp V(pT N -t{U k ) + (1 - V)(pT N -t{u k+1 ) 

B(yj,l k+ i) \ 

+ (j) t (vipT N -t(u k ) + (1 - v)lpT N -t{uk+l)) 

+ (ipt{vipT N -t(uk) + (1 - v)ipT N -t{uk+i)),x)\ (6.23) 

Note that the moment generating function is again exponentially-affine in 
the initial value Xq = x. Here, the moment generating function in (|6,23p is 
well defined for all v G J' k , where 

J k = {veR: vi> t {u k ) + {l-v)i> t {u k+l )eT T , t€ [0,T k ]}. 
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Concluding, we have that all forward prices are of exponential-affine form 
under any forward measure and the model structure is always preserved; as 
a consequence, the model is analytically tractable in the sense of Axiom 2 
with respect to all foward measures. 

Remark 6.3. Note that for the model to make sense and be easy to use 
and implement we must know the functions 4> an d tp explicitly, and not only 
as e.g. the (unknown) solution of a Riccati ODE. 

Remark 6.4. Note that the dynamics of forward prices under their forward 
measure are particularly simple, avoiding recursive formulas, due to the 
specification of the martingales directly with respect to ~Pt n ■ 



7. Interest rate derivatives 

The most liquid interest rate derivatives are caps, floors and swaptions; 
in practice, LIBOR models are typically calibrated to the implied volatility 
surface of caps and at-the- money swaptions, and then hedging strategies and 
prices of exotic options are derived. Thus, it is important to have "closed 
form" valuation formulas for caps and swaptions, so that the model can be 
calibrated in real time. Here we derive such formulas for caps and swaptions, 
making use of Fourier transform methods. 

Caps are series of call options on the successive LIBOR rates, termed 
caplets, while floors are series of put options on LIBOR rates, termed floor- 
lets. Caplets and floorlets are usually settled in arrears, i.e. the caplet with 
maturity T k is settled at time T k+ \ = T k + 5; for simplicity we consider a 
tenor structure with constant tenor length 5, although this assumption can 
be easily relaxed. A cap has the payoff 

N-l 



J2$(L(T k ,T k )-K)+. (7.1) 



fc=i 



Keeping the basic relationship (|2.ip in mind, we can re-write caplets as call 
options on forward prices: 

5(L(T k ,T k ) - K) + = (1 + 5L(T k ,T k ) - 1 + 5K)+ 



Mrp 



X) , (7.2) 



where % := 1 + 5K. 

Each individual caplet is typically priced under its corresponding forward 
measure to avoid the evaluation of a joint law or characteristic function; in 
our modeling framework we have that 



C(T k , K) = B(0, T k+ i) Ep Tfe+i [S(L(T k ,T k ) - K) 

'fe+i 



£(0,T fe+1 )E PT _ (jj^gr-K) 



k 



Then, we can apply Fourier methods to calculate the price of this caplet as 
an ordinary call option on the forward price. 
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Proposition 7.1. The price of a caplet with strike K maturing at time T k 
is given by the formula 



u , ^A k +B k -x Tk {R ~ iv) ^ 



{R-iv){R- 1-iv) 



(7.4) 



where R £ J k n (l,oo) and the moment generating function h.A k +B k -x Tk is 
given by (|6.23p via 

W*^ (») = ^|^^ + , [e'^ B ^] . (7.5) 
Proof. Starting from (|7.3p and recalling the notation (|6.2ip . we get that 

^A k +B k -X T 



C(T k ,K) = B(0,T k+1 )TE^ 



k+1 



X 



(7.6) 



hence we can view this as a call option on the random variable A k + B k ■ Xx k ■ 
Now, since the moment generating function AA k +B k x Tk of A k + B k ■ Xx k is 
finite for R E J' k , and the dampened payoff function of the call option is 
continuous, bounded and integrable, and has an integrable Fourier transform 



for R G (l,oo), we can apply Theorem 2.2 in Eberlein et al. (2009) and 
immediately get that 



C(T k ,K) = gt°^*±l) j 



E P \ e (R-™)(A k +B k x Tk )i 

T k + 1 L J 

(iv - R)(l + iv - R) 



-dv, 



which using (|6.23p yields the required formula. 



□ 



Now, we turn our attention to swaptions, but we restrict ourselves to one- 
dimensional affine processes as driving motions. Recall that a payer (resp. 
receiver) swaption can be viewed as a put (resp. call) option on a coupon 
bond with exercise price 1; cf. section 16.2.3 and 16.3.2 in Musiela and 
Rutkowski (119971) . Consider a payer swaption with strike rate K, where the 
underlying swap starts at time Tj and matures at T m (i < m < N). The 
time-Tj value is 

§ Ti (K,T h T m ) = (l- jr c k B(Ti,T k )\ , (7.7) 

V k=i+i J 

where 

f K, i + l<k<m-l, , . 

Ck = \ 1 + K, k = m. (7 ' 8) 

Now, we can express bond prices in terms of the martingales M u , as follows: 
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since the product is again telescoping. Analogously to forward prices, cf. 
(|6.7p . the dynamics of such quotients is again exponentially affine: 



■j^r = exp [4)T N -t{u k ) ~ (f>T N -t{Ui) 

+ (^T N -t{Uk) - 1pT N -t(ui),X t ] 

=: exp (A kji + B kii ■ X t ). 



(7.10) 



Then, the time-0 value of the swaption is obtained by taking the Prr 
expectation of its time-Tj value, hence 



§ (K,T u T m ) = B(0,Ti)V PTi 
= 5(0,T,)Ep Ti 
= J B(0,T i )Ep T . 



1- Y, CkB(Ti,T k ) 

>, k=i+l / 



Ck 



k=i+l 



1 ~ E c k e Ak ^ +Bk ^ XT * 

^ k=i+l 



, (7.11) 



and this expectation can computed with Fourier transform methods. 
Define the function / via 



f{x) = 1 " E e fc e^e^-* . (7.12) 

V k=i+i J 



We will also assume that, at least some, initial LIBOR rates are positive. 

Proposition 7.2. The price of a swaption with strike rate K , option ma- 
turity T{ and swap maturity T m is given by 



§ (K, T it T m ) = B ^ 2 J i ' I A XTi (R - iv)f(v + iR)dv, 



where the Fourier transform of the payoff function f is 



S(v + iR) = e 



I m a A k i+B k i~y 
(iv-R)y [ k 



/-^ B ki -R + iv iv - R 
\k=i+i ' / 



(7.14) 



Here y denotes the unique zero of the function f , the W^-moment generat- 
ing function Ax T . of X^ is given by (|6.15p and R £ Z l n (0, oo). 
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Proof. Starting from (|7.1ip and using Theorem 2.2 in Eberlein et al. (2009) 
again, we have that 



S {K,Ti,T m ) = B(0,Ti)E Fl 



1 - ^2 c k e Ak > i e Bk < rXT i 



k=i+l 



f(l- jt c k e A ^e B ^A P Ti ,x Ti (dx) 

p V k=i+l / 

Ax T XR-iv)f(v + iR)dv, (7.15) 



5(0, Ti) 
2ir 



where Ax T . denotes the -moment generating function of the random 

variable X^, and / denotes the Fourier transform of the function /. 

Now, we just have to calculate the Fourier transform of / and show that 
the prerequisites of the aforementioned theorem are satisfied; we know that 
Ax^ is finite for R £ Z\ 

Regarding the Fourier transform of /, define the function /, where 

m 

f(x) = 1 - E CfceWw*; 

k=i+l 

since we assumed that some LIBOR rates are positive, Proposition 16.11 and 
Lemma 14.21 yield that B k ^ > for some k. Hence, we can easily deduce 
that f'(x) > 0, therefore / is a strictly increasing function. Moreover, it 
is continuous and takes positive and negative values, hence it has a unique 
zero, which we denote by y. Therefore, 



f{x) =l(x)l { y t00 y 
Now, for z £ C with > 0, the Fourier transform of / is 

f( z ) = J e izx (l - jr c k e Ak -e Bk -A dx 

o \ k=i+l / 



(7.16) 



y 

oo 

/ 



' m \ 

1 - Yl c k e Ak ''e Bk >'' x I dx 

v k=i+l ) 



y 



e izx dx - c ^ Ak ' 1 \ e {iz+Bk ^ )x dx 

y 



k=i+l 



E 

\k=i+l 



c k e 



A k>i +B k>i y 



B k A + iz 



1 

iz 



(7.17) 



Moreover, by examining the weak derivative of the dampened payoff function 
g{x) = e~ Rx f(x) for R > 0, we see that it is square integrable, as is g itself. 
Hence g lies in the Sobolev space i7 1 (lR) and applying Lemma 2.5 in Eberlein 
et al. (|2009p yields that the Fourier transform of g is integrable. □ 
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Remark 7.3. The above valuation formula can be re-expressed as 

§ (K, Ti,T m ) = jB( ° 2 J Af) J Ax Ti (R-iv)f(v + iR)dv, (7.18) 

R 

where 

A XTi (R-iv) = ^l^- A XTi (R-iv). (7.19) 
8. Examples 

Here we present four concrete specifications of the affine LIBOR model 
we have constructed; the driving processes are a Cox-Ingersoll-Ross pro- 
cess and three OU-type processes driven by Levy subordinators, namely the 
Gamma subordinator, the inverse Gaussian subordinator, and a compound 
Poisson subordinator with exponential jumps, such that the OU process has 
the Gamma law as stationary distribution (cf. Nicolato and Venardos 120031 
IKeller-Ressel and Steiner 2008p . We first describe or construct the driving 
affine processes with values in M^o an d then discuss the affine martingales 
used to model LIBOR rates. In the case of the CIR driving process we 
derive a closed- form pricing formula for caps and swaptions, using the x 2 - 
distribution function. 

8.1. CIR martingales. The first example is the Cox-Ingersoll-Ross (CIR) 
process, given by 

dX t = -\{X t -6)dt + 2r 1 ^/x t dW u Io = i£% (8.1) 

where X,9,rj £ K^o- This process is an affine process on M^Oi with 

F(u) = \9u and R(u) = 2r?V - \u. (8.2) 

Its moment generating function is given by 

E, [e uXt ] = exp + x ■ ^(«)) , (8.3) 

where 

Mu) = -^log{l-2 V 2 b(t)u) and ^(«) = ^-^^L- (8.4) 



with 

, . . It, if A = 

b(t) = < ' - xt , and a(t) = e~ xt 



The martingales defined in (J5JJ) thus take the form 

Mf = exp (<fr N - t (u) + (ip TN -t(u),X t )) 

6XP [ ~ V bg (1 " 27]2b(TN " t)u) + l-Wb{T N -t)u • Xt 

.5) 

where u must be chosen such that u < 2rpb{T N ) • ^°^ e * na ^ 7a = oo (see 
Definition I5.5p . such that by Proposition 16.11 the model can fit any term 
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structure of initial LIBOR rates. 



In order to describe the marginal distribution of this process, we de- 
rive some useful results on an extension of the non-central chi-square dis- 
tribution; we say that a random variable Y has location-scale extended 
non-central chi-square distribution with parameters (fi,a,u,a), or short 
Y ~ LSNC— x 2 (a*> °~, v i <y)i 

if ^ has non-central chi-square distribution 
with parameters v, a. The density and distribution function of Y can be de- 
rived in the obvious way from the density and distribution of the non-central 
chi-square distribution. We will also need the cumulant generating function 
of Y, which is given by 

v . . aau 1 

For any •d < we may consider the random variable Y$ with distri- 
bution function F$, defined through the exponential change of measure 
<Uj£ _ g a;i?-«;(i9) _ j^. j g we jj known that the cumulant generating function 

of Y# is given by k$(u) = k(u + i?) — For the LSNC— x 2 -distribution a 
simple calculation using (|8.6|) shows that 



Y$ ~ LSNC-x 2 (ft, ^, v, ^) , with C = 1 - 2<ti? > 0. (8.7) 



Let us now return to the CIR process X. Comparing (|8.4p and (|8.6p . 
shows that 

X, ^ LSNC-x 2 (o, ^6(t), ^.^), (8-8) 

i.e. the marginals of X have LSNC— % 2 -distribution under the terminal mea- 
sure. By (|6.16|) ~ (|6,18|) we know that the measure change from the terminal 
measure to the forward measure Pr fe+1 is an exponential change of 

measure, with & = ipT N -t(u k+ i). Thus, we derive from (|8.7|) that 

X PTfc+1 T^NP v 2 fn X ° Xa{t) \ (RQ\ 

V C(t,T k+1 ) rf rj 2 b(t)((t,T k+ i) J 

where 

C(t, Tfc+i) = 1 - 27? 2 6(t)^ TiV -t(n fc+ i). (8.10) 
Finally, it follows from ()6.2ip . that the log- forward rates have distribution 

log (* {t > Tk \) LSNC-x 2 (A k , B ^ Xe Xa ® 



B(t,T k+1 )J V C(t,r fc+ i)' T? 2 ' r,%(t)at,T k+1 ) 

(8.11) 

under the corresponding forward measure, where A k and B k are given by 
(|6.2ip . Hence, log-forward rates are LSNC— x 2 -distributed under any for- 
ward measure with different parameters a and a, due to the different £. 

We are now in the position to derive a closed-form caplet valuation for- 
mula for the CIR model. Denoting by M = log ( b^^t^) ) ^h e log- forward 
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rate, it holds that 

C(T kt K) = £?(0,T fc+1 )Ep Tfc+i [(e M - X) + ] (8.12) 

= B(0,T k+1 ) {Ep Tfc+i [e A/ l {M > log3C} ] -XTP Tk+1 [M > logX}} 

= B(0, T k ) W Tk [M > log %}-% B(0, T k+l ) P Tfc+1 [M > log X] , 

where we have used (|6,11|) and X = 1 + 5K. The probability terms can be 
evaluated through the distribution function of the LSNC— ^-distribution. 
After some calculations, we arrive at the following result: 

'logX-A k 



£(T k ,K) = B(0,T k ).xt, ai 



logX-A k \ 2 

X.V,Ct2 



a-2 



(8.13) 

where 30* := XB(0,T k+1 ) and xl, a i x ) ■= 1 ~ xl, a ( x ), with xl, a (x) the 
non-central chi-square distribution function, 

X9 B k r, 2 b{T k ) xa(T k 

VT Cl,2 

and 



«1,2 



ri 2 b(T k )Ci, 2 ' 



(i = ((T k ,T k ), ( 2 = C(T k ,T k+ i). 

In a similar way, a closed-form pricing formula for swaptions can be de- 
rived; by (|7.1ip . using (|6.14p . we get 

S (K,Ti,T m ) 



B(0,Ti)TE FTi 

b(o,t n )-e Fti 

B(0,T n ){e f . 



1 



2^ Cfc M^ 



k=i+l 



T, 



s, fc=i+l 



M Tthx Tl >y } 



c ^ 

k=i+l 



M T- l {x Ti >y}\ f, 



(8.14) 

where y is defined as in (|7.16p . Using the known distribution function of 
X?i under Pj; cf. (j8.9j) . the exponential change of measure formula ([57 
and ()6.14p once again, we arrive at 



So(K, Ti,T m ) = B(0, I}) ■ xl, ai (^M) " E CkB ^ Tfc ) ' OV**) 

fc=i+l 



.15) 



where 



— , cr fc = — ^ and a k 



xa(Ti) 



Remark 8.1. Notice that models based on the one-dimensional CIR process 
are complete in their own filtration due to the continuity of paths and the 
Markov property. The hedging strategy is given by A-hedging. 
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Remark 8.2. A particular feature of models based on the one-dimensional 
CIR process is that the LIBOR rate L(t,Tk,Tk + i) is bounded from below 
by 4(exp^4fc — 1). This is undesirable, but a negligible failure, since usually 
the quantity is close to 0. 

8.2. T-OU martingales. The second example is an OU-process on Mj>o 
such that the limit law is the Gamma distribution. Consider the SDE 

dX t = -XX t dt + dH t , Xq = x £ M^o, (8.16) 

where A > 0. The driving Levy process H = (H t )t>o is a compound Poisson 
subordinator with cumulant generating function 

kcp (u) = - J ^, (u < a) 8.17 
a — u 

where a, j3 > 0; hence, H is a compound Poisson process with jump intensity 
\/3 and exponentially distributed jumps with parameter a. The moment 
generating function of H is well defined for u G Zcp = (— oo,a). The limit 
law of this OU process is the Gamma distribution T(a,/3), i.e. it has the 
cumulant generating function 

Kr ( u ) = _ / 31n(l__) ; (8.18) 

v a 



cf. Theorem 3.15 in Keller-Ressel and Steiner (2008) We call the resulting 
affine process the T-OU process. 

The moment generating function of the random variable Xt, using Lemma 
17.1 in |Sato (1999)1 is 

1E X [e uXt ] = exp K C p{e- Xs u)ds + e~ xt u ■ x) ; (8.19) 
now, using the change of variables y = e Xs and f , 1 , , > da; = — 4 In I ax+b \ 

' ° ° a J x(ax+b) b I x I ' 

we get: 

K CP (e- Xs u)ds = / y—ds = f3 In , 8.20 

o Jo a - e~ As u \ a — u J 

since u € (— oo,a). Hence, the moment generating function in (|8.19|) is 

E x [e uXt ] = exp (p In f a ~ e ~ Xtu \ + e -xt u . ^ j (8 _ 21) 



a — u 

which yields that X is an affine process on D = M^o with 

cj) t Cu) = /3 In (?—^ -\ and ipAu) = e~ xt u, (8.22) 

\ a — u J 

and the functions F and R have the form 

Flu) = and = -An. (8.23) 

a — u 



Therefore, the affine martingales constructed in (|5.ip take now the form 

= exp (<fr N - t {u) + (ipT N -t(u), X t )) 



exp [ ^ ( g ~ e- A(r ""^ | + e -A(i*-t) u . Xf | ? (g 24) 



a — u 
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where u must be chosen such that u £ Iqp H M^o = [0,a). Moreover, we 
have that 7x = oo, hence the model can fit any term structure of initial 
LIBOR rates. 

8.3. OU-r and OU-IG martingales. The third example is driven by 
the Gamma subordinator; let H = (Ht)t>o be a Gamma process, i.e. a 
Levy subordinator where H\ ~ T(a,f3), a, (3 > 0. The cumulant generating 
function of the Gamma variable is 

ii 

kt(u) = -/31n 1 - -), (8.25) 

v a 

and is well defined for u < a. The resulting affine process will be called the 
OU-r process. 

The fourth example is driven by the inverse Gaussian (IG) subordinator; 
let H = (Ht)t>o be an inverse Gaussian process, i.e. a Levy subordinator 
where H\ ~ IG(<5, 7), 5,7 > 0. The cumulant generating function of the 
inverse Gaussian variable is 



(u) =8j- - 2u, (8.26) 



2 , 



and is well defined for u G (— 00, \\. The resulting affine process will be 
called the OU-IG process. 

In the sequel, since the construction of the corresponding affine processes 
and the affine martingales is common for both driving subordinators, we 
simply refer to them as the subordinator and subordinated OU process. 

The subordinated-OU process is an affine process X = (Xt)t>o with state 
space D = M^Oj that satisfies the SDE 

dX t = -X(X t - 0)dt + dH t , X = xe 1R >0 , (8.27) 

where A, 9 > 0. The conditional moment generating function of the subordi- 
nated-OU process is given by 

E x [e uXt ] = exp (o(l - e~ xt )u + J K(e" As u)ds + x • e" A *^ . (8.28) 

Hence, we immediately get that 

(j) t (u) = 9(l-e- xt )u+ f K{e- Xs u)ds (8.29) 

Jo 

and 

^(u) = e~ xt u, (8.30) 

which naturally yield the functions F and R: 

F{u) = \6u + k(u) and R(u) = -An. (8.31) 

Moreover, the set X in which the cumulant generating function is well defined 
is respectively 

7 2 

2p = (— 00, a] and Xig = (— 00, — ]. (8.32) 
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Figure 1. Implied volatility surface for the CIR martingales. 



Therefore, the affine martingales constructed in (|5.ip take now the form 
= exp (<h-t{u) + (ifr- t (u),X t )) 



T-t 



= exp ^(l-e- A(T " i) )n+ J K{e~ Xs u)ds + e'^-^u ■ X t j , (8.33) 

where u must be chosen such that u £ I D M^o- 

9. Numerical illustration 

In order to showcase some prototypical volatility surfaces resulting from 
the proposed models we consider the tenor structure of zero coupon bond 



prices from the Euro zone on February 19, 2002; cf. Table [T] and Kluge 
(|2005t pp. 50). We fit the initial LIBOR rates implied by the bond prices 
using the u's as described in Proposition 16. 11 and then price caplets and plot 
the implied volatility surfaces for different parameters of the driving affine 
factor process. The strikes we consider range from 1% to 6% with step 0.5%. 
The implied volatility surface corresponding to the CIR parameters 

A = 0.001, 9 = 0.50, rj = 0.59, X = 1.25 

is shown in Figure [TJ It is important to point out that the closed-form 
caplet formula for the CIR model is several times faster than the Fourier- 
based pricing formula for the same model; on the same computer system 
and using similar implementation the closed form method takes less than ^ 
sec, compared to more than 10 sees for the Fourier method. 

In the example for the r-OU process we consider the same tenor structure 
and the strikes range from 2.5% to 7% with step 0.5%. The implied volatility 
surface corresponding to the r-OU parameters 

A = 0.01, a = 2.00, = 1.00, X = 1.25 
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T 


0.5Y 


1Y 


1.5Y 


2Y 


2.5Y 


B(0,T) 


0.9833630 


0.9647388 


0.9435826 


0.9228903 


0.9006922 


T 


3Y 


3.5Y 


4Y 


4.5Y 


5Y 


B(0,T) 


0.8790279 


0.8568412 


0.8352144 


0.8133497 


0.7920573 



Table 1. Euro zero coupon bond prices on February 19, 2002. 




Figure 2. Implied volatility surface for the T-OU martingales. 

is shown in Figure Let us note that the Fourier pricing formula works 
faster for the T-OU model than for the CIR model, yielding results in about 
7 sees for the whole surface. 

Appendix A. Construction of Markov martingales 

Let X = [Xt)t>o be a time-inhomogeneous Markov process on a general 
state space D, e.g. D = M. d x IR| , starting from x G D; let P Sj t denote its 
transition function from s to t. Consider a "good" function, e.g. / G Cb(D), 
and define the function 

F(t,T,x) = -E x [f(X T )\T t ). (A.l) 

Then, the process (F(t, T, X t )) 0<t<T is a martingale. Indeed, apart from 
finiteness, is suffices to show the martingale property; using the transition 
function and applying the Chapman-Kolmogorov equation, we have 

m x [F(t,T,X T )\T s ] = P s>t P t , T f{X s ) = P s , T f(X s ) 
= F(s,T,X s ). 

In case X is a time-homogeneous Markov process, then we can define the 
function 



F(t,x) = E x [f(X t )}, 



(A.2) 
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and the process (F(T — t, X t )) 0<t<T is a martingale. 



References 

Andersen, L. and R. Brotherton-Ratcliffe (2005). Extended LIBOR mar- 
ket models with stochastic volatility. J. Comput. Finance 9, 1-40. 

Black, F. (1976). The pricing of commodity contracts. J. Financ. Econ. 3, 
167-179. 

Brace, A., D. Gatarek, and M. Musiela (1997). The market model of 

interest rate dynamics. Math. Finance 7, 127-155. 
Duffie, D., D. Filipovic, and W. Schachermayer (2003). Affine processes 

and applications in finance. Ann. Appl. Probab. 13, 984-1053. 
Eberlein, E., K. Glau, and A. Papapantoleon (2009). Analysis of Fourier 

transform valuation formulas and applications. Appl. Math. Finance. 

(forthcoming, arXiv/0809 . 3405). 
Eberlein, E. and F. Ozkan (2005). The Levy LIBOR model. Finance 

Stock. 9, 327-348. 

Filipovic, D. (2005). Time-inhomogeneous affine processes. Stochastic Pro- 
cess. Appl. 115, 639-659. 

Glasserman, P. and S. G. Kou (2003). The term structure of simple for- 
ward rates with jump risk. Math. Finance 13, 383-410. 

Jacod, J. and A. N. Shiryaev (2003). Limit Theorems for Stochastic Pro- 
cesses (2nd ed.). Springer. 

Jamshidian, F. (1997). LIBOR and swap market models and measures. 
Finance Stoch. 1, 293-330. 

Jamshidian, F. (1999). LIBOR market model with semimartingales. Work- 
ing Paper, NetAnalytic Ltd. 

Kallsen, J. and A. N. Shiryaev (2002). The cumulant process and Esscher's 
change of measure. Finance Stoch. 6, 397-428. 

Kawazu, K. and S. Watanabe (1971). Branching processes with immigra- 
tion and related limit theorems. Theor. Probab. Appl. 16, 36-54. 

Keller- Ressel, M. (2008). Affine processes - theory and applications to 
finance. Ph. D. thesis, TU Vienna. 

Keller- Ressel, M. and T. Steiner (2008). Yield curve shapes and the as- 
ymptotic short rate distribution in affine one- factor models. Finance 
Stoch. 12, 149-172. 

Kluge, W. (2005). Time-inhomogeneous Levy processes in interest rate 
and credit risk models. Ph. D. thesis, Univ. Freiburg. 

Miltersen, K. R., K. Sandmann, and D. Sondermann (1997). Closed form 
solutions for term structure derivatives with log-normal interest rates. 
J. Finance 52, 409-430. 

Musiela, M. and M. Rutkowski (1997). Martingale Methods in Financial 
Modelling. Springer. 

Nicolato, E. and E. Venardos (2003). Option pricing in stochastic volatility 
models of the Ornstein-Uhlenbeck type. Math. Finance 13, 445-466. 

Papapantoleon, A. and M. Siopacha (2009). Strong Taylor approxima- 
tion of SDEs and application to the Levy LIBOR model. Preprint, 
arXiv/0906.5581. 



A NEW APPROACH TO LIBOR MODELING 



29 



Revuz, D. and M. Yor (1999). Continuous Martingales and Brownian 
Motion (3rd ed.). Springer. 

Sandmann, K., D. Sondermann, and K. R. Miltersen (1995). Closed form 
term structure derivatives in a Heath-Jarrow-Morton model with log- 
normal annually compounded interest rates. In Proceedings of the Sev- 
enth Annual European Futures Research Symposium Bonn, pp. 145- 
165. Chicago Board of Trade. 

Sato, K. (1999). Levy Processes and Infinitely Divisible Distributions. 
Cambridge University Press. 

Siopacha, M. and J. Teichmann (2010). Weak and strong Taylor meth- 
ods for numerical solutions of stochastic differential equations. Quant. 
Finance, (forthcoming, arXiv/0704 . 0745). 

Walter, W. (1996). Ordinary Differential Equations. Springer. 

Department of Mathematical Methods in Economics, Vienna University of 
Technology, Wiedner Hauptstrasse 8-10, A-1040 Wien, Austria 
E-mail address: {mkeller .papapan, jteichma}@f am. tuwien. ac . at 



